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Abstract 

The OPERA collaboration [Ij has reported evidence of superluminal z/^ propagation between 
CERN and the LNGS. Cohen and Glashow [2j argued that such neutrinos should lose energy 
by producing photons and e^e^ pairs, through TP mediated processes analogous to Cherenkov 
radiation. In terms of the parameter 5 = iv'^—v'^)/v^, the OPERA result corresponds to (5 ~ 5-10~^. 
For this valu^of 6 a very significant deformation of the neutrino energy spectrum and an abundant 
production of photons and e"'"e~ pairs should be observed at LNGS. We present an analysis based 
on the 2010 and part of the 2011 data sets from the ICARUS experiment, located at Gran Sasso 
National Laboratory and using the same neutrino beam from CERN. We find that the rates and 
deposited energy distributions of neutrino events in ICARUS agree with the expectations for an 
unperturbed spectrum of the CERN neutrino beam. Our results therefore refute a superluminal 
interpretation of the OPERA result according to the Cohen and Glashow prediction [2] for a weak 
current analog to Cherenkov radiation. In particular no superluminal Cherenkov like e"'"e^pair or 
7 emission event has been directly observed inside the fiducial volume of the bubble chamber like 
ICARUS TPC-LAr detector, setting the much stricter limit of 6 < 2.5 • 10~* at the 90% confidence 
level, comparable with the one due to the observations from the SN1987a 



Note that {v^ - v^)lv^ « | « 2.5 • lO"'^ 
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I. INTRODUCTION 



The OPERA collaboration has presented evidence of superluminal neutrino propaga- 
tion [T] , reporting a travel time between CERN and the LNGS laboratory some 60 ns shorter 
than expected for travel at light speed. The OPERA result corresponds to 5 = {.'^'t~'^1) hi ~ 
5 ■ 10~^ with only small variations over the energy domain of the detected neutrinos. Obser- 
vations of neutrinos from Supernova SN1987a at much lower energies around 10 MeV yield a 
strong constraint ^ 5 < 4 ■ 10~^ implying a rapid increase with energy of the hereby alleged 
anomaly. 

As is well known, charged particles travelling at speeds exceeding that of light emit char- 
acteristic electromagnetic radiation known as Cherenkov radiation. Because neutrinos are 
electrically neutral, conventional Cherenkov radiation of superluminal neutrinos does not 
arise or is otherwise weakened. However neutrinos do carry electroweak charge and, as 
pointed out by Cohen and Glashow [2], may emit Cherenkov-like radiation via weak inter- 
actions when traveling at superluminal speeds. Cohen and Glashow argue that, under the 
assumptions of the usual linear conservation of energy and momentum and only slow vari- 
ation of 5 over the OPERA-relevant energy domain, superluminal neutrinos would radiate 
and lose energy via the three following processes 

I'x I'x+l (1) 
I'x + l^y + l^y (2) 

I'x ^ i^x + + e~ (3) 

The emission rate and energy loss is dominated by the third process, which is kinemat- 
ically allowed under the stated assumptions. The process [3} from now on referred to as 
pair bremsstrahlung pj, proceeds through the neutral current weak interaction and has a 
threshold energy Et^r ~ 2me/ a/5 corresponding to about 140 MeV for the OPERA value of 
5. In the high energy limit the electron and neutrino masses may be neglected, and Cohen 
and Glashow [2J comput^the rate of pair emission F, and the associated neutrino energy 
loss rate dE/dx to leading order in 5: 

35 1927r3 " ^ ' 

dx" 112 1927r3^'^^ 



^ These expressions have corrected a numerical factor error in ^ of 4/5. 



Note that the average fractional energy loss per pair emission event is ^^^^ ~ 0.78; that is, 
about I of the neutrino energy is lost on average with each emission. Furthermore, under 
the approximation of a continuous energy loss, the integration of dE/dx over a distance L 
provides the following result for the final neutrino energy, E^f, as a function of the initial 
energy, E^i. 

^ ^ - ^ (6) 



Elf El 112 1927r3 

Folding the initial neutrino spectrum of the CERN to Gran Sasso neutrino beam with the 
energy at Gran Sasso predicted with the above formula, the expected neutrino interaction 
rates and pair bremsstrahlung rates as a function of 5 may be estimated. In particular, for 
5 = 5- 10^^ and L = 732 km, equation |6] would predict that few neutrinos with energy larger 
than ^13 GeV would reach Gran Sasso. 

However, since neutrinos lose a large fraction of their energy at each pair creation event, 
and the resulting deflection angles are not negligible with respect to the angular width of 
the CNGS neutrino beam, a continuous energy loss approximation is suitable only for a 
qualitative estimate of the spectral distortion. 

A full three-dimensional Monte Carlo calculation of the propagation of neutrinos from 
CERN to Gran Sasso has therefore been performed for several values of S (see section [ll|). 
As a result, the expected rates of neutrino charged current interactions and of e^e~ pair 
events have been obtained as a function of 6 and are presented in the next section. 

The rates obtained in this way have been compared with the results of year 2010 and 
part of the year 2011 exposures of the reconstructed neutrino charged current events in the 
ICARUS/CNGS2 experiment located in Hall B of the Gran Sasso Laboratory. 



II. SIMULATION RESULTS 

A full 3-dimensional simulation of the generation and transport of CNGS neutrinos from 
CERN to Gran Sasso while undergoing pair bremsstrahlung has been performed using the 
official CNGS simulation setup [HIITO], based on the fluka fU] 112] Monte Carlo transport 
code. 

Accounting for the threshold, and under the hypothesis that S does not vary significantly 
in the range of energies of interest, the pair bremsstrahlung interaction rate differential in 
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TABLE I. Expected neutrino and e^e rates at Gran Sasso. All rates are given for a 1 kt detector 
and 10^^ pot. 





cc 


NC 


CC>60 GeV 




s 


(all flavours) 


(all flavours) 


(i^M + i^/i) 







644 


203 


57 





5 • 10-8 


644 


203 


57 


27 


5 • 10-7 


643 


203 


56 


2.1 • 10^ 


5 • 10-6 


594 


188 


8.5 


7.2 • 10^ 


5 • 10-5 


203 


85 


< 10-6 


1.1 • 10^ 



the neutrino energy loss, w, and in the pair invariant mass Sg+e- = can be expressed as: 
d^r G%5^ s so\i 



dwds 1927r3 El V s 
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The kinematical limits are given by 



(7) 

So = ^ (8) 



E;>w^>s>So (9) 



The neutrino deflection angle ^ with respect to the incident neutrino direction can be 
expressed as: 

cos^^l- ^'-f + (10) 



and the e+e pair angle as: 



The resulting mean free path for a 19 GeV neutrino (the fluence-averaged energy of CNGS 
neutrinos) is ~ 490 km for 5 = 5- 10-^, and the deflection angle is of the order of com- 
parable with the angular width of the neutrino beam. Hence the need for a full Monte Carlo 
simulation of the neutrino propagation to Gran Sasso. 

All results presented in this section are for 10^^ protons on target (pot) and, for rates, for 
a detector (Argon) mass of 1 kt. In this way they can be easily rescaled to whichever Gran 
Sasso detector mass and exposure, neglecting the minor differences in neutrino cross sections 
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FIG. 1. Computed spectra at Gran Sasso for 6 = (solid line, green), and for (5 = 5 • 10 
(dashed line, blue). The units are u cm~'^ GeV^^ lO^^^pot^^ , where pot=protons on target. 



in an energy range dominated by DIS among Argon and other materials. The nominal yearly 
number of protons for CNGS is 4.5 ■ 10^^ pot. The statistical error on integrated values (e.g. 
total rates, total fluence, etc) is less than one percent in all cases. The systematic error on 
the computed neutrino (and hence e"'"e~ pairs) rates is mostly due to the uncertainties in 
the hadron production model of FLUKA, and can be conservatively estimated to be lower 
than 10% (see for example pTl fT3] ). 

The unperturbed {6 = 0) fluence spectra of CNGS t'^ at Gran Sasso, and the one com- 
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FIG. 2. Computed e^"e^ pair spectra at Gran Sasso for 5 = 5- 10"^ (solid line, green), 5 = 1- 10~^ 
(dashed line, blue, multiplied by 1000), and 6 = 5- 10~^ (dotted line, red, multiplied by 10^). The 
event rate units are GeV~^ for a 1 kt detector and 10^^ protons on target (pot). 

puted corresponding to S = 5 ■ 10~^ are shown in Fig. [Tj the lack in the latter spectrum of 
the sharp 12.5 GeV ridge predicted by formula [6] can be easily appreciated. 

The computed spectra of the expected events due to e"^e~ pairs at Gran Sasso are shown 
in Fig. m for (5 = 5 ■ 10-^ 1 ■ 10"^ 5 ■ 10"* respectively. 

The computed (anti)neutrino charged and neutral current rates (all flavours included), the 
charged current rates for i/^ and z/^ with energy above 60 GeV, and the pair bremsstrahlung 
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FIG. 3. Computed neutrino charged current rates (dashed hne, red, left vertical scale), and 
Vfj, charged current rates (x 100), > 60 GeV (dotted line, green, left vertical scale), e+e" rates 
(solid line, blue, right vertical scale) at Gran Sasso, for 10^^ protons on target (pot) and 1 kt 
detector. 



rates at Gran Sasso are presented in Fig. [3j The expected rates are summarized in Table [T] 
for a few representative values of 5. 
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III. EXPERIMENTAL RESULTS AND RELATED CONSTRAINTS 




6.2cm, 66.8cm = ni ± 26 MeV/c 

= 28.0 ± 2.5" 

FIG. 4. Typical event recorded in ICARUS. Evidence for a pair of 7's from a vr" (tracks 16a and 
16b) with a momentum of 912 MeV/c pointing at the primary vertex, showing the typical behavior 
of 7 conversions in the TPC-LAr Imaging chamber. 

The ICARUS experiment [5H8] consists of 760 t of super-pure Liquid Argon operated as 
a very high resolution Time Projection Chamber, namely a Bubble Chamber-like detector 
recording all events with an energy deposition in excess of a few hundred MeV within a 
window of 60 /xs centered around the neutrino pulse from the CERN-SPS. The hadronic and 
electromagnetic energy depositions of each event are accurately measured by calorimetric 
determination while the muon momenta are measured with the help of the multiple scattering 
along the very many points of the long muon tracks. An example of an event with a pair of 
7's produced by a secondary 7r° is shown in Fig. |4| 

In the following an analysis of neutrino interactions from the 2010 and part of the 2011 
CNGS runs is compared with the expectations for different values of 6. While the fine details 
of the analysis could still somewhat evolve, it will be shown that 6 values of the order of the 
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<EData> =10.1 ± 0.5 GeV 
rms =9.6 GeV 



< Ej^c > =10.1 ± 0.2 GeV 



rms =11.6 GeV 



< Ej^c 5=5.10 '' > = 5.2 ± 0.05 GeV 
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FIG. 5. Experimental raw energy E^ep distribution for and i/^ CC interactions in ICARUS 
(blue symbols) compared with the Monte Carlo expectations for an unperturbed spectrum (red 
solid histogram), and for 6 = 5- 10~^ (green dashed histogram). 



one claimed by OPERA can be readily excluded on the basis of the observed rates and raw 
energy deposition spectra. In order to carry out as much as possible a bias-free analysis, the 
raw energy deposition distributions as recorded in ICARUS with the calorimetric methods 
have been compared with the expectations of a full Monte Carlo simulations of the detector: 
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FIG. 6. Experimental (blue symbols) raw energy deposition distribution for neutral current events 
compared with the Monte Carlo expectations (red histogram) for the unperturbed CNGS spectrum. 
Only experimental and Monte Carlo events with energy deposition greater than 500 MeV have been 
considered. 



only the correction for the signal quenching in LAr has been applied to both the experimental 
and Monte Carlo results. 

A dedicated search for e+e~ events has been also carried out using the same exposure. 
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TABLE II. Observed and expected neutrino and e'''e~ rates at Gran Sasso for the ICARUS ex- 
periment. Both the experimental and computed rates are normahzed to the exposure used for the 
analysis of the corresponding experimental channel (see text for details). 



Rates 


Observed 


Expected 


Expected 






6 = 


6 = 5- 10-5 


CC 


308 


315±5 


98.1±2 


NC 


89 


93.1±3 


33.0±1 


CC, Edep > 25 GeV 


25 


18 ±1.3 


< 10-6 


e"'"e~ 








7.4 • 10*^ 



This analysis constrains 6 to values a few order of magnitude smaller than the one claimed 
by OPERA. 

The number of collected neutrino interactions has been compared with the predictions for 
the CERN SPS neutrino beam in the whole energy range, corrected for the fiducial volume 
and DAQ deadtime. The experimental analysis corresponds to an integrated exposure of 
6.70 ■ 10^^ t-pot. This exposure is the combination of a fiducial volume of 447 t of Liquid 
Argon and a total number of protons on target (pot) at CERN of 4.9 ■ 10^^ for the year 2010, 
and 434 t and 1.04 ■ 10^^ pot for the fraction of the year 2011 analyzed up to now. This 
exposure applies fully to the search for pair bremsstrahlung events which do not require any 
further cut on the fiducial volume because of their expected clean signature. 

In order to identify z/^ and z/^ charged current (CC) as well as neutral curent (NC) events, 
further cuts have been applied to the fiducial volume: in particular events with the vertex 
in the last 2.5 metres of the detector have not been considered for this analysis in order to 
identify cleanly possible muon tracks. The total number of identified z/^ and CC events, 
and of neutral current events, are compared in Table [Tij 21 events cannot be safely assigned 
despite the reduced fiducial volume. The resulting reference exposure for NC and CC events 
is 5.05 ■ 10^1 t-pot. 

The measured raw energy deposition E^ep for and CC events as obtained from a 
calorimetric measurements corrected only for signal quenching is presented in Fig. |5] The 
experimental distribution is compared with a full Monte Carlo simulation of the experimental 
apparatus for 5 = and 6 = 5-10-5. The experimental finding matches well the Monte Carlo 
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expectations for the unaffected CERN neutrino beam. 

Tfie analysis of neutral current events gives similar results. The experimental and ex- 
pected spectra of the deposited energy for NC events are presented in Fig. |6j only events with 
experimental or Monte Carlo energy deposition in excess of 500 MeV have been considered, 
in order to avoid possible misidentifications or inefficiencies. 

The strong constraints of Cohen and Glashow [2] predict that a superluminal high energy 
neutrino spectrum will be heavily depleted and distorted after L = 732 km from CERN to 
LNGS: in particular for 6 in the range indicated by OPERA the charged current rate would 
be reduced to roughly 32% of the expected one, the average z/^ energy would be 12.1 GeV 
(against 19 GeV), the average energy of undergoing charged current interactions would be 
12.5 GeV (against 28.7 GeV), and no neutrino interactions should be observed above 30 GeV. 
Indeed at 5 = 5 ■ 10^ essentially no E^, > 30 GeV should arrive from CERN to LNGS, while 
our results are indicating no visible deviation of the incoming neutrino beam with respect to 
the expected rate and energy distribution. This result confirms the inconsistence between 
the OPERA 6 value and the observed neutrino rate and spectrum already reported in ref . [2] . 

In addition, with ICARUS a bubble chamber like detector a much more stringent limit to 
6 may be set from the direct observation inside the ICARUS detector volume of Cherenkov 
generated by the passing superluminal neutrinos. These events would be 
characterized either by a single gamma ray converting into an e^e~ pair (eq. [T]) and/or two 
single electrons (eq. |3]) both with no hadronic recoils in the incoming neutrino direction. The 
transverse momenta of the particles in the events ([T]) and ([3]), as determined by the centre 
of mass system, are however far too small to be experimentally observable. Therefore events 
of both types ([!]) and ([s]) would appear as narrow e"*" e~ pairs pointing directly to the beam 
direction, with no detectable hadronic activity. The rate of such events for the ICARUS 
detector exposure under consideration (6.70- 10^^ t-pot) can be derived from those presented 
in Fig. [3] With the OPERA result (5 ~ 5 • 10^^) more than 7 ■ 10^ electron positron pairs 
should have been observed for this exposure, each with an energy spectrum peaked around 
10 GeV (see Fig. [2]). This number should not be a surprise since the CC neutrino event rates 
of the 2010 and the 2011 ICARUS samples with average energy of ~28 GeV represents a 
total of 2.1-10^^ incoming neutrinos. No Cherenkov like event has been detected in ICARUS. 
The experimental event rates are compared in table |TT] with the expected ones, for 6 = and 
6 = 5- 10^^. The results presented in the table clearly show that the latter value for 6 can be 
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like events (eq. Tp 



excluded. Taking into account both the absence of narrow e^e~ pairs pointing directly to the 
beam direction and the presence of several high energy charged current events (for instance, 
25±5 events with deposited energy in excess of 25 GeV are present in Fig. [5] against 18±1.3 
expected), we can set the limit 5 < 2.5 ■ 10"^ at 90% CL for CNGS neutrino^ comparable 
to the limit 6 < 1.4 ■ 10~® established by SuperKamiokande |i3j from the lack of depletion 
of atmospheric neutrinos, and somewhat larger than the lower energy velocity constraint 
6 < A - 10~^ from SN1987a [1]. The ICARUS events already collected during 2011 represent 
conservatively a factor three higher statistics and should provide more accurate information 
on the indicated process. A similar increase in statistics is expected from the 2012 CNGS 
run. However, due to the 5^ dependence of the pair bremsstrahlung cross section, no major 
change of the 6 limit can be expected if no e+e^ event will be found in the final data sample. 

CONCLUSIONS 

The spectra and rates at Gran Sasso of neutrino and e"'"e~ for the CNGS beam have 
been computed in the theoretical framework presented in [21 [T3]. In particular, pair 
bremsstrahlung events have been accounted for during the propagation of neutrinos from 
CERN to Gran Sasso National Laboratory. The resulting neutrino spectra and rates for 
6^5- 10~^ as suggested by OPERA are significantly different from the unaffected ones. 
Preliminary results from the ICARUS experiment do not support any statistically significant 
deviation from the unperturbed spectrum and therefore exclude 6 values comparable to the 
one claimed by OPERA. 

Furthermore ICARUS did not detect any e+e" event, despite a few millions were expected 
for (5 = 5 • 10"^ The lack of e+e"-like event translate into a 90% CL limit of 5 < 2.5 ■ 10~^ 
for multi-GeV neutrinos. 
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